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ABSTRACT 

Atmospheric  density  fluctuations  induce  phase 
perturbations  that  degrade  the  spatial  coherence  of  a  laser 
beam  propagating  through  the  atmosphere.  These 
degradations  spread  the  laser  beam  and  alter  the  centroid 
and  intensity  profile  stochastically.  Turbulent  conditions 
are  found  at  virtually  all  levels  of  the  atmosphere.  A 
substantial  fraction  of  the  optical  turbulence  along  a 
vertical  path  arises  from  the  heat  flux  between  the 
atmosphere  and  the  Earth's  surface.  This  type  of 
turbulence  is  typically  within  the  first  100  to  200  meters 
above  the  surface. 

During  this  thesis  research,  a  high  frequency  acoustic 
sounder  was  developed  to  analyze  this  turbulent  layer.  The 
primary  focus  was  the  development  of  the  command  and 
control  software  required  to  coordinate  the  data  collection 
and  reduction.  The  system  was  used  at  two  sites  and  should 
prove  useful  in  quantifying  the  effects  of  optical 
turbulence  within  the  surface  boundary  layer  on  laser  and 
optical  system  performance. 
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I.   INTRODUCTION 

A  coherent  laser  beam  propagating  through  the 
atmosphere  is  very  susceptible  to  numerous  turbulence 
inflicted  degradations.  As  electromagnetic  waves  transit 
the  turbulent  regions,  atmospheric  irregularities  randomize 
the  amplitude  and  phase  of  the  wave.  In  order  to  quantify 
the  altitude  dependence  of  the  atmospheric  turbulence,  a 
high  resolution  turbulence  profiler  is  needed.  Acoustic 
echosounders  are  frequently  used  to  detect  and  measure 
atmospheric  density  and  velocity  irregularities  resulting 
from  air  currents,  temperature  inversions,  humidity 
variations,  mechanical  turbulence  and  other  density 
fluctuations . 

Presently  two  atmospheric  optical  parameters ,  the 
spatial  coherence  length  ( ro )  and  the  isoplanatic  angle 
(9o)  are  measures  of  the  perturbation  of  an  electromagnetic 
wave  propagating  through  the  atmosphere,  and  are  accurately 
measured  by  optical  systems  developed  by  Walters  [Refs.  1 
and  2] .  Although  these  systems  measure  a  path  integral  of 
atmospheric  turbulence  with  high  accuracy,  a  major  drawback 
of  each  system  is  that  no  provision  is  available  for 
determining  the  height  of  the  atmospheric  disruptions.  If 
these  disturbances  are  found  to  exist  very  near  the 
surface,   it  may   be  possible   to   negate   their   effect  by 


elevating  the  optical  systems  or  controlling  the  generation 
of  turbulence. 

This  thesis  deals  with  the  design,  construction  and 
implementation  of  a  high  frequency,  acoustic  echosounder 
which  will  accurately  analyze  the  atmospheric  density 
fluctuations  within  approximately  200  meters  of  the 
surface.  As  this  project  is  a  product  of  the  research  and 
efforts  of  two  students,  the  work  was  appropriately 
divided.  My  particular  task  was  to  devise  a  computer 
program  which  would  control  all  aspects  of  the  echosounder 
operation.  Areas  of  particular  interest  were  the 
integration  of  hardware  and  software  within  the  data 
acquistion  system,  developing  plotting  algorithms, 
controlling  the  system  timing,  and  setting  the  input 
parameters  which  determine  the  range  and  sensitivity  of  the 
device.  My  colleage,  LT .  Weingartner  [Ref.  3],  dedicated 
his  efforts  toward  the  actual  design  and  hardware 
development  of  the  echosounder. 

Acoustic  echosounders  have  been  developed  and  in  use 
for  many  years  and  have  proven  to  be  valuable  probes  for 
analyzing  the  structure  and  dynamics  of  the  lower 
atmosphere  [Refs.  4  through  7].  Devices  similar  to  ours 
have  been  used  to  obtain  profiles  of  the  atmospheric 
density  and  temperature  fluctuations  [Refs.  8  and  9].  Our 
device  used  a  high  speed  HP217  computer  to  control  and 
monitor    the    echosounder    which    provided    real   time 
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information.  Additionally,  we  had  the  ability  to  store  and 
reproduce  the  atmospheric  profile  plots  at  will.  This 
information  was  crucial  in  analyzing  the  atmospheric 
turbulence  degradations  on  laser  and  electro-optical 
system's  performance. 

This  document  shall  address  the  theoretical  background 
of  acoustic  sounder  operation  in  Chapter  II.  Chapter  III 
shall  present  a  synopsis  of  the  system  design  and 
associated  software.  A  summary  of  actual  data  collected  at 
two  sites  is  analyzed  in  Chapter  IV  and  Chapter  V  discusses 
the  conclusions  and  recommendations  found  as  a  result  of 
our  research. 


II.   BACKGROUND 

Acoustic  echosounders  probe  the  atmosphere  by 
transmitting  a  pulse  of  acoustic  power  which  is 
subsequently  scattered  back  from  the  atmosphere  by 
temperature  and  velocity  inhomogeneities .  The  echosounder 
(echosonde)  equation,  often  referred  to  as  the  radar 
equation  in  meteorology,  is  used  to  determine  the 
backscattered  acoustic  power.  This  equation  is  summarized 
by  Neff  in  Reference  10  and  is  based  upon  the  work  of 
Tatarski  [Ref.  11]  and  Little  [Ref.  5]. 

Pr    =    Er  [    PtEt     ][    exp(-2aR)    ][    ao(R.f)    ][    ^cx    ][    AGR-  2     ] 

where 

Pr  is  the  electrical  power  returned  from  a  range  R, 
Pt  is  the  electrical  power  transmitted  at  frequency  f , 
Er  is   the  efficiency  of  conversion  from  acoustic  power 

to  electrical  power  by  the  transducer, 

Et  is   the   efficiency   of   conversion   from  electrical 

power  to  acoustic  power  by  the  transducer, 

exp(-2aR)   is    the   round    trip   power   loss   due   to 

attenuation  where  a  is   the  average   attenuation  ( meters- i ) 

to  the  scattering  volume  at  the  range  R  (meters), 

ao(R,f)  is  the  scattering  cross  section  per  unit  volume 
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at  a  distance  R  and  frequency  f , 

c  is  the  local  speed  of  sound  (meters/second), 

x    is  the  pulse  length  (seconds), 

A  is  the  aperture  area  of  the  antenna  (meters2), 

R  is  the  range  (meters),  and 

G  is  the  effective-aperture  factor  of  the  antenna. 

Empirically  measuring  or  calculating  the  values  for  all 
other  terms,  one  can  use  this  equation  to  determine 
cro(R,f),  the  scattering  cross  section  per  unit  volume;  that 
is,  the  fraction  of  incident  power  backscattered  per  unit 
distance  into  a  unit  solid  angle  at  a  frequency  f.  Based 
upon  experimental  results,  Tatarski  [Ref .  11]  expresses  the 
acoustic  backscatter  cross  section  per  unit  volume, 
ao(R,f),  in  the  equation, 


Ct2 

cro(R,f)  =  0.0039  ki/3   

To2 


where 

k   =  2k/  A    is   the   incident   acoustic   wavenumber  at 

wavelength  A, 

To  is  the  local  mean  temperature  in  degrees  Kelvin,  and 
Ct2  is  the  temperature  structure  parameter. 

Combining   this   equation   with  the  echosonde  equation,  one 

obtains  a  volume-averaged  measure  of  Ct2. 
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11  2  1  PR 

CT2    = To2    k-i/3 R2exp(2aR) 

0.0039      ErEt  ct        AG        Pt 


Hall  and  Wescott  [Ref.  12]  calculated  a  beam-shape 
compensation  factor  of  0.40  for  a  piston  source  antenna 
with  a  uniformly  illuminated  square  aperture.  This  value 
is  the  same  as  the  effective-aperture  factor,  G,  and  can  be 
substituted  into  the  above  equation.  Approximating  the 
aperture  area  of  the  antenna  to  be  equivalent  to  25  times 
the  aperture  area  of  a  single  speaker  having  a  diameter  of 
7.620  centimeters,  we  get  a  value  of  0.1140  square  meters. 
Combining  these  values  with  the  numerical  constants  in  the 
above  equation,  we  can  simplify  the  equation  for  Ct2. 


1  1  Pr 

Ct2     =    11245 To2    k-l/3 R2exp(2aR) 

Er  Et  ct        Pt 


The   acousto-electrical    efficiency   factors   for   our 
echosounder  were  measured  with  a  calibrated   microphone  and 
an  anechoic   chamber.    The  calculated  values  for  Er  and  Et 
are  each  0.5.   Substituting   these   values   into   the  above 
equation,  we  can  further  simplify  the  equation  for  Ct2  . 


To2       PR 

Ct2    =    44980    k-i/3 R2    exp(2aR) 

ex        Pt 
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The  above  equation  was  incorporated  into  the  computer 
program  that  provides  the  system  control,  data  acquistion 
and  data  processing  techniques  [Appendix  A] .  The-  reduced 
data  was  then  used  to  measure  the  temperature  structure 
parameter  as  a  function  of  time  of  day  and  altitude  for 
various  sites. 
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III.   SYSTEM  DESIGN  AND  EQUIPMENT  DEVELOPMENT 

A.   HARDWARE 

Most  of   the  hardware   used  in  this  system  was  standard 
scientific   equipment   and  is  illustrated  in  Figure  1.   The 
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pieces  which  were  specially  designed  by  LT  Weingartner 
I Ref .  3]  are  the  acoustic  array  and  the  enclosure  for  this 
array.  The  design  considerations  for  the  speakers'  and  the 
array  format  are  given  in  Appendix  B  and  the  design  of  the 
enclosure  to  house  this  array  is  outlined  in  Appendix  C. 
Aside  from  these  two  pieces  of  equipment,  there  are 
primarily  six  other  components  which  comprise  this 
echosounder.  Below  is  a  brief  description  of  each  of  these 
additional  components. 

1- •   HP200  Series  Computer 

The    Hewlett-Packard    (HP)   200   Series   Computer 
includes  a  20  megabyte  hard  disk  drive,  a  floppy  disk  drive 
and  an   associated  printer   and  monitor.    The  HP200  Series 
Computer  was  the  central   control  component   for  the  entire 
echosounder  arrangement.    The   computer  used   was  an  HP::: 
programmed  in  Basic  3.0  and  equipped  with  an   Infotek  BC!  ' 
Basic   Compiler    and   an    Infotek   FP210   Floating   Point 
A  :ceierator  to  enhance  the  speed  of  program  execution. 
program   defined   all   parameters   for  the  HP3314A  Function 
Generator  as  well   as   execute   the   trigger   command  which 
■  c  ' iced  the   echosounder  transmitted   acoustic  pulse.   The 
computer  also  received  the  data  from  the  acoustic  array  via 
the   pre-amplif ier ,   bandpass   filter  and  analog  to  digital 
converter.   The  computer  then  conducted   the  data  reducti 
routines    to    produce    and   display   a   high   resoiuti  n 
atmospheric  profile. 


2.  HP3314A  Function  Generator 

The  Hewlett-Packard  (HP)  3314A  Function  Generator 
was  a  multimode  function  generator  capable  of  providing 
sine,  square  and  triangular  wave  functions  as  well  as  any 
desired  waveform  ranging  in  frequency  from  0.001  Hertz  to 
19.999  Megahertz.  The  HP3314A  Function  Generator  was  used 
to  supply  a  pulse  of  an  integer  number  (usually  100)  of 
sinusoidal  cycles  of  constant  amplitude  to  the  QSC  Model 
1700  Audio  Amplifier.  A  constant  frequency  setting  of  5000 
Hertz  was  used  for  all  data  runs. 

3.  QSC  Model  1700  Audio  Amplifier 

The  QSC  Model  1700  Audio  Amplifier  is  a  high  power 
amplifier  which  can  supply  350  watts  over  an  8  ohm  load. 
This  amplifier  was  used  to  boost  the  output  voltage  of  the 
function  generator  by  a  factor  of  20  from  1 . 5  volts  to  30 
volts.  The  power  supplied  to  each  speaker  in  the  array 
during  operation  was  then  37.5  watts. 
4 •   Pre-amplif ier 

The  pre-amplif ier  was  designed  and  constructed  by 
Walters  [Ref .  13]  to  supply  a  gain  of  1000  to  the  returned 
signal  and  isolate  the  data  acquistion  components  from  the 
transmitted  pulse.  The  pre-amplif ier  may  be  thought  of  as 
a  safety  and  switching  mechanism  for  the  system.  An  LT 
1037  and  an  LT  1007  Operational  Amplifier  were  selected  for 
use  in  the  pre-amplif ier  based  upon  their  low  noise 
properties  which  were  evaluated  by  physically  incorporating 
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them  into  the  device  and  measuring  their  noise  and  gain 
characteristics  with  an  HP3561  Spectrum  Analyzer.  The 
diodes  used  are  IN  4000  Series  Rectifiers  capable  of 
carrying  one  ampere.  These  rectifiers  serve  to  isolate  the 
power  amplifier  from  the  amplifier  input  stage  and  to  limit 
the  voltage  applied  to  the  operational  amplifiers.  The 
electrical  diagram  of  this  device  is  outlined  Figure  2. 


TO  ARRAY 
o 


iokn 


Fig.  2.  Schematic  of  Pre-amplif ier 
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5.   Rockland  Wavetek  Model  852  Filter 

The  Rockland  Wavetek  Model  852  Filter  operated  as  a 
48  db  per  octave  bandpass  filter  to  suppress  the  broadband 
noise  of  the  system.  High  and  low  bandpass  settings  of 
5500  Hertz  and  4500  Hertz  respectively  were  used.  The 
filter  response  at  these  settings  is  illustrated  in  Figure 
3. 
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Fig.  3.  Bandpass  Filter  Response 

The  three  db  (half  power)  bandwidth  is  approximately  1  KHz 
wide.  Although  it  would  be  desirable  to  reduce  the 
bandpass  to  between  50  and  100  Hertz  without  impairing 
signal  response,   the  Rockland  Wavetek  Model  852  Filter  was 
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the   best   filter   available.     It  performed  adequately  in 
spite  of  the  large  lKHz  bandwidth. 

6 .   Infotek  AD2QQ  Analog  to  Digital  Converter 

An  Infotek  AD200  12  bit  Analog  to  Digital  Converter 
was  used  to  digitize  the  signal  voltage  for  the  computer 
system  at  a  12.5  to  20  KHz  sample  rate. 

B .   SOFTWARE 

The  software  and  HP200  Series  Computer  system  were 
responsible  for  controlling  and  monitoring  every 
operational  phase  of  the  hardware  components.  The  program 
which  accomplished  this  task  was  entitled  "ACRDR"  and  is 
listed  in  Appendix  A.  This  program  was  modeled  after  a 
program  written  by  Walters  [Ref.  1.3]  but  each  program 
performs  a  distinctly  different  computational  task.  This 
program  was  written  in  HP  Basic  3.0  and  was  compiled  by  an 
Infotek  BC203  Basic  Compiler  to  enhance  the  speed  of 
execution. 

The  program  "ACRDR"  is  easily  broken  down  into  a  number 
of  blocks  and  subroutines  which  performed  specific 
operations.  These  sections  are  outlined  in  a  flowchart 
(Fig.  4).  Each  block  is  straightforward  in  its  purpose, 
and  the  program  is  designed  to  be  as  helpful  to  the  user  as 
possible.  As  a  prologue  to  the  actual  program  code,  there 
is  a  listing  of  all  the  program  variables  with  a  short 
description  of  their  use.  Such  a  listing  familiarizes  a 
user   with   the   computations   to  be  made  and  also  provides 
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quick   references 
implemented. 


for   any   future   modifications   to   be 
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Offset 

Range 

Pulse  Length 

Attenuation  Coefficient 


I 


Prepare  Screen  for 
Data  Display 


I 


Trigger  Pulse  from 
HP  331 4A 


Send  Cj  Plot  to  the 
Printer 


Collect  Return  and 
Perform  Reduction 
Routines,  Store  Data  if 
Required 


Compute  Values  of  C-^ 
and  Plot 


Send  Plot  to  the 
Printer 


Compute  the  Average 
Noise  Level 


Fig.  4.  Flowchart  of  Fundamental  Program  Operations 

The  initialization  procedures  which  set  all  the 
parameters  necessary  for  data  collection  follow  the 
variable  definition  section.  These  parameters  control  such 
things  as  the  contrast  between  background  and  return  signal 
and  the  setting  of   the   computer's   internal   clock.    All 
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options  of  the  HP3314A  Function  Generator  are  decided  by 
the  user  at  this  phase.  Another  initialization  involves 
preparing  the  Infotek  AD200  Analog  to  Digital  Converter  for 
operation.  Finally,  internal  arrays  are  dimensioned, 
integer  variables  are  defined  and  the  computer  function 
keys  are  redefined  to  suit  the  "ACRDR"  program. 

Now  that  the  system  is  ready  to  operate,  certain 
quantities  must  be  determined  based  upon  atmospheric 
conditions  at  the  time  of  data  collection.  From  the 
ambient  temperature,  a  speed  of  sound  in  air  is  determined 
by  the  equation  in  Kinsler,  et  al ,  [Ref .  14'.  p.  106] 

C  =  Co  [  1  +  (T/273)  ]l/2 , 

where 

C  is  the  local  speed  of  sound, 

Co  is  the  speed  of  sound  in  air  at  0°  Celsius  (331.6 
m/s ) ,  and 

T  is  the  local  temperature  in  degrees  Celsius. 
Using  this  speed  with  a  user  decided  sampling  frequency  and 
number  of   data  points,   the  maximum  range  of  the  traces  is 
determined  by  the  relation, 

RANGE  =  C/2  •  l/( SAMPLE  FREQUENCY)  •  N  , 

where 

RANGE  is  the  range  of  the  echosounder  in  meters, 
SAMPLE  FREQUENCY  is  the  user  desired  rate  of  data 
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acquisition,  and 

N  is  the  number  of  data  points  per  trace. 
Typically,  the  values  of  many  of  the  parameters  are  not 
varied.  The  value  of  N  above  is  almost  always  taken  to  be 
the  allowed  maximum  of  16301  data  points  and  the  speed  of 
sound  can  be  safely  estimated  to  be  340  m/s  at  room 
temperature.  By  using  a  typical  sample  frequency  of  20000 
Hz,  a  range  of  approximately  135  meters  is  obtained. 

Another  parameter  determined  is  the  atmospheric 
attenuation  coefficient.  This  calculation  is  done  in  a 
subroutine  obtained  from  Reference  15  and  is  necessary  in 
this  program  for  the  calculation  of  Ct  2  .  Other 
computations  made  prior  to  program  execution  include  the 
pulse  length,  wavenumber,  and  the  D-C  offset  of  the 
equipment . 

At  this  point,  the  computer  is  finally  ready  for  data 
collection.  The  screen  setup  displays  a  distance  versus 
time  plot  and  the  internal  clock  of  the  computer  is 
synchronized  with  the  pulse  from  the  HP3314A  Function 
Generator.  Following  transmission  of  the  pulse,  the  return 
signal  is  received,  digitized,  stored  in  an  array  and  data 
reduction  commences.  A  block  averaging  technique  is  used 
in  which  a  block  of  data  points  the  size  of  the  number  of 
emitted  cycles  is  summed  and  averaged.  Returns  are  plotted 
as  darkened  areas  on  the  echosounder  traces  with  the 
intensity   of   the   darkened   area   linearly  related  to  the 
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magnitude   of   the   return   signal.   This  same  procedure  is 
used  in  the  computation  of   Ct 2   with   some   very  important 
differences.      Aside   from    simply   correcting   for  nr* 
spherical  divergence,  Ct 2  is  corrected  for  electronic  gain, 
the   ratio   of   power   returned   to   power  transmitted,  the 
efficiency  of  the  speakers   as  transmitters   and  receivers, 
the  area  of  the  speaker  array,  the  atmospheric  attenuation, 
the  pulse  length,   the   temperature,   the   scattering  cross 
section  per   unit  volume   at  a  specific  range  and  frequency 
and  finally  the  effective   aperture  factor   of  the  antenna. 
Additionally,   Ct2 (Range)   is   averaged   for   a   particular 
altitude  over  15  minute   intervals.     Finally,   after  each 
pulse  is   reduced,  a   corresponding  mean  square  noise  level 
is  determined.     This   is   done   by   averaging   each  block 
average   at   maximum   range   until  at  least  ten  values  have 
been  used  in  the  average.   An  upper  limit  of   five  over  the 
average    noise    figure    (this    corresponds   to   voltage 
fluctuations  on  the  order  of   10-7   volts)   is   set   on  the 
routine   to   avoid   averaging   any  strong  return  signals  or 
anomalies  such  as  passing  aircraft.    After  the   ten  values 
are  averaged,   every  subsequent   pulse  is  averaged  into  all 
the  preceding  noise  levels  and  removed  from  each  subsequent 
return  signal. 

Ultimately,  after  each  15  minute  interval,  or  at  the 
users  request,  the  display  terminal  image  is  printed.  At 
this  time  the  Ct 2  computations  are  conducted  and  plotted  on 
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the  screen  and  printer  as  range  versus  logi o  of  Ct  2  .  Upon 
completion  of  the  printing  of  these  plots,  the  system 
begins  the  data  processing  for  the  next  15  minute  interval. 
There  are  certain  options  built  into  the  program  to 
allow  a  user  to  change  various  aspects  of  operation.  The 
function  keys  allow  the  user  to  change  the  local 
temperature,  sample  frequency  or  intensity  factor  during 
program  execution.  Additionally,  the  user  can  quit  or 
restart  the  program,  print  the  partial  trace  on  the  screen 
or  elect  to  save  a  trace  on  a  floppy  disc.  The  save 
routine  is  invoked  for  the  subsequent  15  minute  interval 
after  the  appropriate  function  key  is  depressed.  Saving  a 
future  trace  may  seem  awkward,  especially  if  a  user  would 
like  to  keep  an  interesting  trace  which  is  presently  on  the 
terminal.  This  problem  cannot  be  readily  solved  unless 
each  trace  is  recorded  to  disc  without  user  intervention. 
At  present,  this  is  not  done  because  only  8  traces  (2  hours 
of  data)  can  be  written  to  a  floppy  disc  before  it  is  full. 
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IV.  Ml 

A.  ECHOSOUNDER  PERFORMANCE 

An  analysis  of  the  acoustic  echosounder  output  was 
conducted  to  determine  the  validity  of  previously 
determined  echosounder  parameters.  The  e- *  decay  time 
constant  for  applied  voltage  was  calculated  in  Appendix  B 
to  be  approximately  900  us  or  for  convenience,  1  ms . 
Reviewing  typical  echosounder  traces  and  Ct 2  plots 
[Appendix  D] ,  it  was  evident  that  the  recovery  time  was 
typically  found  to  be  on  the  order  of  33  ms  or  roughly  six 
meters  past  the  end  of  the  transmitted  pulse  length  of  3.4 
meters.  The  recovery  time  is  consistent  with  the  time 
required  for  the  30  volts  on  the  drivers  to  decay  to  the 
microvolt  level.  A  more  detailed  analysis  of  the  hardware 
is  found  in  Reference  3. 

B.  SITE  EVALUATION 

Echosounder  data  was  collected  at  two  different 
locations.  The  primary  data  collection  site  was  the  upper 
roof  of  Spanagel  Hall  at  the  Naval  Postgraduate  School, 
Monterey,  California.  This  site  was  chosen  simply  for 
convenience.  Data  gathered  at  this  location  is  believed  to 
represent  the  California  coast  during  the  spring  near  sea 
level.  The  second  site  chosen  was  in  the  vicinity  of  the 
24   inch    telescope    at    Lick    Observatory,   San   Jose, 
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California.  This  site  is  located  atop  Mt .  Hamilton  at  an 
altitude  of  approximately  5700  feet  and  nearly  20  miles 
inland  from  the  coast. 

These  two  data  collection  sites  represent  areas  of 
differing  atmospheric  air  pressures,  water  vapor  pressures, 
local  temperature  ranges,  and  local  wind  velocity  ranges. 
These  characteristics  all  play  important  roles  in  effecting 
the  local  atmospheric  turbulent  conditions  and  thereby  the 
atmospheric  structure  parameter,  Ct2. 

In  addition  to  collecting  echosounder  data  at  Lick 
Observatory,  simultaneous  measurements  of  the  isoplanatic 
angle  (80)  and  spatial  coherence  length  ( ro )  were  made  with 
systems  developed  by  Walters  [Refs.  1  and  2].  A  basic 
knowledge  of  these  two  systems  is  necessary  to  understand 
the  correlation  procedures  made.  The  isoplanatic  angle 
(80)  is  primarily  an  upper  atmospheric  measurement  which 
indicates  atmospheric  disruptions  at  a  range  of  2  to  15 
kilometers.  The  spatial  coherence  length  ( ro )  is  a  measure 
of  the  effects  of  the  entire  atmospheric  blanket  on 
coherent  light  transmission.  A  close  comparison  of  all 
three  data  sets  should  give  us  an  accurate  description  of 
both  the  lower  and  upper  troposphere  as  well  as  the 
stratosphere  above. 

Based  upon  isoplanatic  angle  (80)  and  spatial  coherence 
length  ( ro )  measurements  at  Mt .  Wilson  in  California,  a 
strong  correlation  between  the   two  measurements   occurs  if 
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the  low  altitude  boundary  layer  contribution  is 
sufficiently  small.  This  strong  correlation  helps  to 
reinforce  the  overall  description  of  the  atmosphere  at  the 
time  of  data  collection.  Figures  5  and  6  graphically 
illustrate  the  atmospheric  measurements  made  at  Mt .  Wilson 
on  2  April  1987  by  Walters. 
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Fig.  5.  Isoplanatic  Angle  Measurements,  Mt .  Wilson 
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Fig.  6.  Spatial  Coherence  Length  Measurements,  Mt.  Wilson 

The  strong  correlation  between  the  isoplanatic  angle 
(9o)  and  the  spatial  coherence  length  (ro )  is  especially 
evident  between  the  hours  of  0700  and  1300  universal  time. 
The  close  tracking  of  these  two  measurements  during  this 
time  interval  indicate  that  the  upper  atmospheric 
conditions,  as  measured  by  9o ,  are  dominating  the  entire 
atmospheric  profile  as  measured  by  ro .  Unfortunately,  the 
existence  or  non-existence  of  any  turbulent  surface  effects 
cannot  be  ascertained  by  the  employment  of  the  above  two 
systems  alone.    However,   employment  of   these  two  systems 
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together  with  the  acoustic  echosounder  should  enable  us  to 
produce  a  complete  atmospheric  profile  with  strong 
correlation  between  all  three  atmospheric  measurements. 

On  9  and  10  April  1987,  all  three  systems  were  operated 
at  Lick  Observatory.  Again,  a  good  correlation  between  the 
isoplanatic  angle  (9o )  and  the  spatial  coherence  length 
( ro )  measurements  was  noted.  In  addition,  a  strong 
correlation  between  the  spatial  coherence  length  ( ro )  and 
echosounder  measurements  was  present.  A  comparison  of  the 
atmospheric  data  in  Figures  7  and  8  shows  a  good 
correlation  between  the  two  parameters  especially  during 
the  0930  to  1130  time  interval  on  9  April  1987.  However, 
during  subsequent  hours  the  isoplanatic  angle  (9o) 
measurements  remain  high  (=12  urad )  indicating  relatively 
calm  turbulent  conditions  in  the  upper  atmosphere  while  the 
spatial  coherence  length  ( ro )  values  drop  sharply  after 
12:00  Universal  Time  indicating  dominant  and  increasing 
lower  atmospheric  turbulence.  This  increase  in  the  low 
level  turbulence  should  be  evident  in  the  echosounder  data 
commencing  around  1200  universal  time  (0400  local  standard 
time)  on  9  April  1987.  A  comparison  of  the  echosounder 
data  in  Figures  9  through  11  illustrates  this  increase  in 
the  local  surface  turbulence. 


29 


9  Rpr I 1  1937  Lick 


Fig.  7.  Isoplanatic  Angle  Measurements,  Lick  Observatory 
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Fig.  8.  Spatial  Coherence  Length  Measurements 
Lick  Observatory 
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Fig.  9.  Echosounder  Trace,  Lick  Observatory 
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Fig.  10.  Echosounder  Trace,  Lick  Observatory 
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Fig.  11.  Echosounder  Trace,  Lick  Observatory 

The  strong  correlation  between  the  echosounder  data  around 
4:00  to  5:00  Standard  Time  and  the  ro  measurements  around 
12:00  Universal  Time  combined  with  the  lack  of  correlation 
between  the  9o  and  ro  measurements  indicate  that  the  lower 
atmospheric  and  surface  turbulence  are  dominating  the 
atmospheric  profile  during  this  time  period. 

Data  collected  on  10  April  1987,  again  illustrate  the 
strong  correlation  between  the  three  atmospheric 
measurements  made.  A  comparison  of  the  9o  and  ro 
measurements  during  the  time  interval  of  0700  and  1300 
Universal  Time  indicate  steady  turbulent  conditions  in  the 
upper  atmosphere   and  greatly   varying  turbulent  conditions 


32 


at  lower  atmospheric  levels.  This  is  evident  in  Figures  12 
and  13  by  the  consistent  values  of  80  during  the  time 
period  compared  with  the  steady  increase  and  eventual 
decline  of  the  ro  values  during  the  same  time  interval. 
The  trace  variations  in  Figure  13  during  the  hours  of  0800 
and  1200  are  indicative  of  a  period  of  decreasing  lower 
atmospheric  or  surface  turbulence  followed  by  the  onset  of 
an  increasingly  turbulent  period  around  1200.  This 
turbulent  trend  is  strongly  supported  by  the  echosounder 
data  in  Figures  14  through  21. 
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Fig.  12.  Isoplanatic  Angle  Measurements,  Lick  Observatory 
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Fig.  13.  Spatial  Coherence  Length  Measurements, 
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Fig.  14.  Echosounder  Trace,  Lick  Observatory 
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Fig.  15.  Echosounder  Trace,  Lick  Observatory 
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Fig.  16.  Echosounder  Trace,  Lick  Observatory 
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Fig.  17.  Echosounder  Trace,  Lick  Observatory 
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Fig.  18.  Echc  ounder  Trace,  Lick  Observaiory 


36 


135 


10   flpr    198? 


LICK   OBSERVRTORY 


en 
u 

u 


u 
c 

-7 

£ 


3:35  3:40 

TIME    LOCRL    (STD) 


3:45 


Fig.  19.  Echosounder  Trace,  Lick  Observatory 
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Fig.  20.  Echosounder  Trace,  Lick  Observatory 


37 


135 


cn 
on 
u 

i- 

Ld 

Z 


u 
o 

-7 

QL 


90 


45 


10   Rpr    1987 

a 

.    :».     ■-. 


LICK  OBSERVATORY 


i.  1-fi 


•'  ::"4  •  •  • 

:  :<  •  •  • 


.  .-.  :.jn-.     .....;*.■  ::yi 


..  *i  ."■:■>■  • 

1  .:•  V  .   *  I 


@^«?w^-ife:H  lit  £   *w.  Mfc  * ; 


^a^ffl'A^wiiiUWMMa 


4:35  4  :  40 

TIME  LOCRL  (STD) 


4:45 


Fig.  21.  Echosounder  Trace,  Lick  Observatory 

The  large  200mm  coherence  lengths  around  11:00  Universal 
Time  are  consistent  with  the  low  turbulence  evident  in  the 
echosounder  profiles  around  3:00  Standard  Time  followed  by 
a  pre-dawn  increase  in  the  surface  turbulence. 

Only  acoustic  sounder  atmospheric  measurements  were 
made  on  the  upper  roof  of  Spanagel  Hall  at  the  Naval 
Postgraduate  School.  Data  runs  on  26  and  27  April  1987, 
included  both  echosounder  data  and  the  associated  Ct 2 
plots.  A  representative  sample  of  the  data  collected 
during  this  period  is  included  in  Appendix  D.  Of 
particular  interest  is  the  data  in  Figure  22,  which  shows 
the  maritime  boundary  inversion  layer  at  about  100m  and 
being  perturbed  by  convective  plumes  at  lower  altitudes. 
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Fig.  22.  Echosounder  Trace  and  Ct 2  Plot,  NPS  Roof 
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V.   CONCLUSIONS  AND  RECOMMENDATIONS 

A  high  frequency  ( 5KHz )  phased  array  echosounder 
constructed  to  measure  low  level  turbulence  appears  to  work 
within  the  10  to  135  meter  altitude  range.  When  using  the 
device,  detailed  profiles  of  the  short  range  atmospheric 
density  fluctuations  were  obtained.  The  profiles  were 
found  to  correlate  very  well  with  the  measurements  of  the 
isoplanatic  angle  (80 )  and  the  spatial  coherence  length 
( ro )  during  periods  of  simultaneous  operation.  This  short 
range  echosounder,  when  used  in  conjunction  with  the  other 
atmospheric  measuring  devices,  is  an  invaluable  tool.  The 
graphical  output  provides  us  with  a  more  complete 
description  of  the  atmosphere  that  can  be  used  to  calculate 
Cn2 ,  the  atmospheric  index  of  refraction  structure 
parameter. 

Further  software  development  such  as  the  incorporation 
of  Fast  Fourier  Transform  routines  into  the  echosounder 
program  will  determine  the  radial  velocity  profile  of  the 
return  signal.  Such  an  improvement  will  determine  the 
velocity  of  the  probed  air  masses  and  plot  these 
echosounder  traces  as  a  function  of  color  intensity. 

Other  areas  of  further  reasearch  include  the  design  and 
testing  of  different  array  patterns.  One  such  pattern,  the 
hexagonal   array,   is   already   undergoing   tests.    As  the 
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speaker  array  evolves,  it  is  also  evident  that  the 
enclosure  must  follow  suit  to  accomodate  the  new  array 
design  and  associated  beam  patterns.  Further  improvements 
in  software  routines  are  also  inevitable.  It  is  always 
desirable  to  store  on  a  disc  all  data  collected  at  a  site. 
Presently,  the  floppy  disc  capacity  is  inadequate  for  the 
storage  of  more  than  two  hours  of  data.  Furthermore,  if 
FFT  and  Doppler  routines  are  to  be  added,  the  computational 
speed  of  the  computer  will  inhibit  the  pulse  repetition 
rate.  A  special  data  acquisition  technique  called  "Direct 
Memory  Access"  may  then  have  to  be  introduced  into  the  HP 
computer  to  allow  simultaneous  data  collection  and 
processing. 

Finally,  this  product  may  be  used  in  various 
applications  not  already  addressed.  One  such  use  may  be  to 
measure  windshear  at  airports.  Other  researchers  have 
expressed  interest  in  using  the  echosounder  to  measure 
arctic  atmospheric  conditions  during  meteorological 
surveys.  Assuredly,  as  the  device  evolves  into  a  more 
compact  and  highly  efficient  instrument,  its  range  of 
application  will  continue  to  expand. 
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APPENDIX  A 


ACOUSTIC  ECHOSOUNDER  PROGRAM 
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***#**#*#*****#»*»#*******#♦*•*»*»*******#*■»****»***##*#**#*## 
»  ACRDR  * 

*  COMPUTER  SOFTWARE  WRITTEN  BT  LT.  M.  WROBLEWSKI  FOR  AN    * 

*  ECHOSOUNDER  BUILT  FOR  A  MASTERS  OF  SCIENCE  THESIS      * 
»       BY  LT.  M.  WROBLEWSKI  AND  LT.  F.  WEINGARTNER         * 

*  ADUISOR:  PROF.  D.  WALTERS  * 

*  JUN.  1387  * 
#******♦♦***#*****#**♦♦*♦»***■»♦»*#♦♦*******#*«#*****«*♦»♦**♦* 


The  computer  program  receives  information  from  an  acoustic 
array  through  an  A-D  converter.   This  information  is  the 
returned  signal  of  an  acoustic  pulse  as  it  passes  through 
the  atmosphere.   The  data  is  then  used  to  display  the  return 
intensity  with  distance  as  a  function  of  time.   This  is  a 
short  range  device  (from  150m  to  200m). 


LIST  OF  VARIABLES 

I  &  K  u  J   -   counters  for  loops 

Rec_num   &   Nrec   -  real  and  integer  representation  of  the 

record  number  for  storage 

Plotnum   -   counter  used  to  insert  form  feeds  between  plots 

Disc auuressj!   ~~   Storage  location  o  i  u3id  fixe 

FileJ 


j  j3 


380 


410 


N I  Sr  Num   ~   string  a n u   i  e a  ±    representation  o  i  the  number  o 
cycles  and  the  number  o i  points  used  in  the 

computation  of  the  uiock.  average 

Point(»)   -   real  number  representing  the  noise  and  range 
corrected  average  over  given  number  of  points 

Dat(-)   -   array  of  a-d  converter  output  after  sampling 
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423 
430 
440 

453 

4G3 

473 

403 

490 

503 

510 

520 

530 

540 

550 

5G3 

570 

500 

530 

G00 

610 

G20 

G30 

G40 

G50 

GG0 

G73 

GO0 

G90 

700 

T  1  n 
I   1  XJ 

720 

I   -it) 

740 
750 
7G0 
770 
700 

-inn 

Oron, 


043 
050 
000 
070 
000 
090 
900 


D2( *  )   -   array  of  reduced  and  averaged  data  including 
offset,  correction  for  noise,  and  range 
corrections 

Hrs   -   the  integer  hours 

Min  -   the  number  of  minutes 

Qtrhr  &  Qtrmin   -   keeps  track  of  the  passage  of  each 

15  minute  interval 

Numl   -   counter  for  computation  of  block  averages 

Num2   -  counter  for  the  number  of  block  averages  made 

M   -   the  counter  used  in  the  plot  label  routines 

Timcnt   -   time  interval  between  data  samples  in 
nanoseconds 

Kmax  -   total  number  of  block  averages  computed 

Freq$  &  Freq   -   the  frequency  input  of  the  HP3314A  in  kilohertz 

Zone!  -   the  appropriate  time  zone  the  operator  desires 

Site!  -   the  name  of  the  appropriate  site  of  data  collection 

Gave_plt  &  Cntrl   -   on/off  toggles  used  to  determine  whether 

a  particular  run  will  be  saved  to  a  disc 

Nplot   -   the  number  of  increment  along  the  entire 
horizontal  axis 

Offset   -   the  computed  D-C  offset  for  the  system  prior 
to  data  acquisition 

Noise   -   a  running  total  of  the  noise  accumulated  at  the 
■  ■  5 ,\ i r u pi  range;  used  i o  rind  and  re m ova  una 
a  v er  a ge  noise 

Limit   -   the  upper  bound  on  the  noise  figure  used  to  insure 
that  large  returns  are  not  included  in  the 
noise  computations 

|<1   -   a  counter  of  the  number  of  traces  used  to  determine 
the  noise  average 
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910 

920 

930 

940 

950 

9G0 

970 

980 

990 

1000 

1010 

1020 

1030 

1040 

1050 

1060 

1070 

1030 

1090 

1  100 

1110 

1120 

1130 

1  140 

1  150 

1  1G0 

1  170 

1  180 

1190 

1200 

1210 

1220 

1230 

1240 

1250 

12G0 

1270 

12G0 

1230 

1303 

i  -7  i  a 
I  j  l<] 

i  7i  n> 
I  _)  l.  a 

1  jj'u 

1340 
1350 
13G0 
1370 
1380 
1390 
1400 


Samave   -   a  sample  average  used  in  the  computation  of  the 
D-C  offset  , 

Off_ave   -   a  toggle  used  to  inhibit  the  firing  of  the 
HP3314A  during  D-C  offset  computations 

Count!  u  Cnt   -   the  string  and  real  representation  of 
the  number  of  data  points  to  be  taken 


!  Gamfreq 

! 

!  Tminc 


the  user  input  sampling  frequency  desired 

1  over  the  sampling  frequency;  the  time  between 
samples 


TimeS   -   the  string  required  by  the  A-D  converter  to 
sample  at  the  desired  rate 

COkelvin   -   the  speed  of  sound  at  0  degrees  celcius  in  m/s 

Temp   -   temperature  of  the  surroundings  in  degrees  celcius 

5pd   -   the  speed  of  sound  in  air  computed  for  the  input 
temperature 

Rspd   -   the  relative  speed  of  sound  that  the  echosounder 
sees  which  is  half  the  computed  speed 

Rdist   -   the  distance  traveled  in  one  time  increment 

Maxrng   -   the  maximum  range  of  the  echosounder  found 

by  multiplying  the  distance  per  time  increment 
by  the  number  of  time  increments 

Exrng   -   the  range  rounded  to  the  nearest  value  evenly 
divided  by  15;  used  solely  for  plotting 
purposes 

Far   -   used  to  estimate  the  far  field;  this  value  keeps  the 
range  correction  from  being  applied  to  data  very  near 
to  the  source 

Ilvl   -   the  U5cr  input  intensity  level  divisor;  this  value 
sets  the  screen  contrast  in  data  return 

3cale_y   -   the  ratio  of  exrng  to  maxrng;  this  value  is  used 
to  keep  the  plotted  range  of  data  accurate 

Npoint  u  Npointl  2<  Npoint2   -   variables  used  to  keep  track  of 

time  passage  between  pulses 
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1410 

1423 

1430 

1440 

1450 

UG0 

1470 

1430 

1490 

1500 

15T0 

1520 

1530 

1540 

1550 

15G0 

1570 

1530 

1530 

1G00 

1510 

1620 

1G30 

1G40 

1G50 

1G50 

1G70 

1  Go0 

1G30 

1700 

1710 

i  nn  rk 
I  I  L.XJ 

i  i7n 
i  i   _  w 

I  I  "»0 

1750 
17G0 


1730 


1  3G0 

I  O  I  <o 
I  oou 

1830 

i  nnrn 

I  -  ■_  - 


T2  S<  T1  St  T0   -  used  to  synchronize  data  collection  with  the 
clock 

TIMEDATE   -   the  internal  clock  of  the  computer 

Dvd   -   divisor  of  the  data  average;  used  because  computer 
multiplication  is  faster  than  division 

R   -   the  range  of  a  block  of  data  samples 

Run_ave   -   the  running  sum  of  the  block  samples 

X   -   the  value  of  the  data  point  less  the  D-C  offset 

Ns   -   the  final  average  of  the  block'  of  data  points;  also 
used  in  the  computation  of  the  noise  figure 

Corr   -   the  noise  correction  applied  to  the  data  samples 

Time   -   the  horizontal  position  of  the  trace  on  the  plot 

Timedist   -   the  horizontal  width  of  the  trace  on  the  plot 

End_pt   -   the  last  point  of  the  plot  vertically  taking 
into  account  the  scaling  factor  scale_y 

Inc   -   The  vertical  increment  along  the  plot 

Z   -   The  final  reduced  data  points  which  are  output  on  the  plot 

Dis   -   The  vertical  height  of  the  trace  on  the  plot 

Ap$ ,Frq$ ,Nm$ ,En$,UoS,Hz$   -   strings  needed  to  set  the  HP3314A 

Amp$   -   userv  input  amplitude  for  the  HP3314A 

Pls_lng   -   the  pulse  length  of  the  burst 

Ct(»)   -   the  atmospheric  temperature  structure  parameter 

Y   ■   the  number  of  traces  u s e u  l n  computing  u  u  % :  ) 

Ypl   -   the  vertical  position  on  the  Ct  plot 

Pnt   -   the  horizontal  position  on  the  Ct  plot 

i 

i 

!  K3   -   the  wave  number  to  the  l/3rd  power 

I 

!  A   -   The  area  of  the  receiver  (array) 
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1310 
1923 
1933 
1943 
1350 
1GE0 
1973 
19S0 
1G33 


L. 


G   -   The  effective  aperature  factor 

Er   -   efficiency  of  coversion  of  acoustical  power 
to  electrical  power  on  the  recieve  side 

Et   -   the  efficiency  of  conversion  of  electrical 

power  to  acoustical  power  on  the  transput  side 

Pt   -   The  computed  transmitted  power  to  the 
acoustic  array- 


Gain   -   the  electronic  gain  of  the  equipment 
Zimp   -   the  speaker  array  impedence 


2033 

2043 
2050 

4.0G  0 

2073 

2380  OPTION  BASE  1 


2133 

2113 

2123 

2133 

c  1  4«j 

2158 

21G3 

2173 

2180 

2190   I 

2200   ! 
>« 

nn-in 


initialize  the  arrays   &   set  dimensions 
declare  all  integer  variables 


DIM  Disc_addressSC20] ,File$C30] .Point ( 300 > ,Ct ( 200 5 
INTEGER  I  ,Hr  .Plotnum  ,Print_key  ,Num .Nurnl  ,Num2  ,M  ,K 
INTEGER  Rec_num .Krnax 
INTEGER  D2(300),Dat( 1S301 )  BUFFER 


initialization  routines 
Generator 


set  time,  set  HP3314A  Function 


2250  Rstrt:CALL  Freq_init < Freq$ ,N$ ) 

22B0  CALL  Imt_ad233 

2273  CALL  Set_t ime( Zone$ ) 

2233  INPUT  "SITE  NAME  ".SlteS 

i.._Ju  Javc_pl  t-</ 

n-rr\n  ri      jl  _  1  _/» 

< "u*j  Uliuli.-*/ 


i.  j  i  <; 

T7-?ni 

i-  -i  U  <J 

n  -7  a  n, 
i it«J 

">  ~*<~  [7\ 

23G0 

TJ7CI 

c  -J  i   <J 

2380 
2390 


keyboai  d  seu up 


s  e  u  s  labels  on  the  computer  function 


!  CLEAR  THE  KEYS 


OUTPUT  KBD;  "SCRATCH  KEYE" ; 

CONTROL  2,2;1 

ON  KEY  1  LABEL  "NEW      TEMP"  GOTO  Speed 
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i.480 
2410 
2420 

1  A  t  f\ 

n  a  4  r* 

L  -  -  - 
ito<! 

2430 

2500 

2510 

2520 

25j0 

2540 

2550 

25G0 

2570 

2550 

2590 

2G00 

2S10 

2S20 

2S30 

2G40 

2G50 

2GG0 

2G70 

2G80 

2G30 

n  —irxn 

L.  t   <J  <J 


-   I 

2743 
2750 
27G0 

L.   I     I  t) 

Tinni 
L.  i  uu 

2730 


ON  KEY  2  LABEL 
ON  KEY  4  LABEL 
ON  KEY  5  LABEL 
ON  KEY  G  LABEL 
ON  KEY  7  LABEL 


' INTENS.    FACTOR"  GOTO  Lvl 
"PRINT    TRACE"  GOTO  Prt_dmp 
'GAVE  NXT  PLOT"  GOTO  Gvpt 
'CHANGE   SMPL  FRQ"  GOTO  Gfrq 
'RESTART"  GOTO  Rstrt 


ON  KEY  0  LABEL  "QUIT"  GOTO  Quit 


!  Set  constants 


Di5c_addres5$=" :  ,700,1  ,0' 

Maxrec=723 

Npiot=300 

Plotnum=1 

Offset=0 

GINIT 

Noise=0 

Linit-10 

Y=0 

K1=0 

Samave=0 

A-.1140 

G=.4 

Er=.49S 

Et=.43G 

Gain=7312G"2 

Zimp»12.12 


Input  and  calculation  of  terms  used  by  the  A-D  converter  and 
the  computation  of  range 


PRINT  "Enter  the  number  of  data  points  desired  (max  1G301)' 
PRINT  "Use  increments  of  the  number  of  cycles  then  add  1" 
PRINT  "ex:  (150  cycles  X  100  data  points/cycle  )+1  =  1 5001" 
PRINT  "RECOMMEND  USING  MAX  VALUE  OF  1G301" 
INPUT  Counts 


A  a  m 


Pt-(Et*(30A2))/<2*Zimp > 


2550  : 

23G0 

2370 

2830 

2330 


r-frq 


PRINT  "Input  the  sampling  frequency  desired 
PRINT  "determine  the  range  of  the  sounder. 

PRINT  "    SAM.  FREQ.  -  12500 RANGE  =  225  M" 

PRINT  "    SAM.  FREQ.  =  20000 RANGE  -  135  M" 

INPUT  "ENTER  THE  SAMPLE  FREQUENCY" .Samfreq 


This  frequency  will 
Examples  are : " 
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2900  GCLEAR 

2910  Tminc=1 ./Samfreq 

2920  Timcnt=l000»INT(Tminc/1 .E-B  ) 

2930  TimeS="TIME  "&<JAL$(  Timcnt  ) 

2940  Cnt=VAL< Counts) 

2950  C0kelvin=331  .S 

29S0  Num-VAL(NS) 

2970  Freq=VAL(FreqS)»1000 

2933  I 

2990  ! 

3020  !   Computation  of  the  D-C  offset  prior  to  program  run 

3010  ! 

3020  ! 

3030  0ff_ave»1 

3040  Offset-0 

3050  FOR  1=1  TO  10 

30G0  CALL  Read_ad200(Dat( *) .Counts, TimeS , Of f_ave) 

3070  FOR  K-l  TO  (Cnt-1 ) 

3080  Samave=5amave+Dat ( K ) 

3090  NEXT  l< 

3100  PRINT  "COMPUTING  D-C  0FF5ET" 

31 10  NEXT  I 

3120  Off_ave=0 

3130  0ffset»Samave/< 10»<Cnt-l  )  > 

3140  PRINT  "OFFSET  IS  :" .Offset 

3150  ! 

31G0  ! 

3170  I   computation  of  the  speed  of  sound  at  a  given  temp  and  the 

3180  !   range  of  detection  of  the  device 

3190  ! 

3200  I 

3210  Speed: INPUT  "Enter  the  temperature  (celsius)  "  ,Temp 

Spd=C0kelvin»(SQR( 1 +< Temp/273 ) ) ) 

Rspd=Spd/2 

Rdist=Tminc*Rspd 

3250  Maxrng=Cnt*Rdi5t 

32G0  Exrng=< ( INT< Maxrng/30 ) )+.5)*30.0 

3270  Far=1 . /( Tminc *Spd ) 

32S0  Lambda=Spd/Freq 

3280  K3-< (2*PI  )/Lambda)"( 1/3) 

3333  Pl5_lng=(R5pd*Num )/Freq 

771    A  PAI     I            AXi_._.._i,_r,/r-___        T        ..  _         /\ll \ 


_  ^  _  , 


Ui.t-1 


[0 


I 


JJJU 

3340  PRINT  "Enter  the  relative  intensity  division  level.   Thi, 

3350  PRINT  "value  is  used  to  determine  the  plot  intensity  by 

33G0  PRINT  "dividing  the  block  average  sum  by  this  number. 

3370  PRINT  "This  value  is  dependent  upon  the  gain  of" 

3380  PRINT  "the  device  and  will  probably  need  adjustment" 

3390  PRINT  "during  run.   Start  with  a  value  of  about  4000" 

3400  PRINT  Pt 
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3410 
3420 
3430 
3440 
j450 
34S0 

-7  *  nrt, 
Jt80 

3430 
3500 
3510 

7r-i  r* 

w'w-  CI 

3530 
3540 
3550 
35S0 
3570 
3580 
3590 
3600 
3610 
3620 
3630 
3640 
3650 
3560 
3670 
3680 
3660 
3700 


3740 
3750 
3760 
3770 
3780 
3790 

-}  nrsrs 
udvo 

TO  1  r\ 


3840 
3850 
3860 
3870 
3880 
3890 
3900 
3910 


Lvl:  INPUT  Ilvl 


!   set  up  the  plot 


Again: CALL  Plot_5etup( Nplot ,Site$ .Maxrng  ,Soale_y ,Exrng ,2one$ ) 
Nrec=0 
FOR  1=1  TO  200 

Ct(I)=0 
NEXT  I 
Y=0 

Npoint 1 =0 
OUTPUT  KBD;"L"! 
! 
! 

Sync:  !   synchronize  data  collection  with  clock 
T1=INT(TIMEDATE  MOD  86400) 
IF  TKT0  THEN  T0=T0-86400 
IF  T1-T0<1  THEN  GOTO  5ync 
T0=T1 


f 


Data  collection  and  reduction 


Read_sig:    !   read  the  A-D  converter 

CALL  Read_ad200(Dat(» ) .Counts ,Tine$ , Of f_ave ) 

Npoint=INT(T1  MOD  3600  MOD  Nplot) 

Npoint 2=Npoint -Npoint 1 

Numl =Nun-1 

Num2=Cnt-Num 

Dvd=l ./(Nun) 
K-0 
! 
! 
FOR  1=1  TO  Num2  5TEP  Nun 
R=I*Rdist 
IF  KFar  THEN  R=1 
!<=K  +  1 
Run_ave=0 
FOR  J=I  TO  I iNunl 

X=Dat( J  )-0f f set 

Run_ave=Run_ave+X  *X 
NEXT  J 

Ns=SQR(Run_ave*Dvd  ) 
Point(K)=R*<Ns-Corr> 
IF  Kl>10  THEN 

Ct(K  )=Ct(K  )+( (Point <K  )"2  )»EXP( 2*R*Atten  )  ) 
Ct(K)  =  (Ct(K)*(  (3.339E-8T2  )*Er  )/Zimp 
END  IF 
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3930 

3340 

3950 

33G0 

3970 

3980 

3990 

4000 

4010 

4020 

4030 

4040 

4050 

40G0 

4370 

4080 

4090 

4100 

4110 

4120 

4130 

4140 

4150 

41G0 

4170 

4180 

4190 

4200 

4210 

4220 

4230 

4240 

4250 

42B0 

4270 

4280 

4290 

4333 

4310 

4323 

4330 

4340 

4350 

43G0 


IF  Point(K  )>ABS(327G7)  THEN  Point  (  K >=327G7 
D2(K  +  10)  =  INT(Point(l<  )  ) 
NEXT  I 
IF  Kl>10  THEN 

Y=Y+1 
END  IF 


Noise  correction  routine 


IF  Ns<(l_imit+5)  THEN 

Noise=Noise+Ns 

Kl-Kl+1 

IF  Kl>10  THEN 

Corr=Noise/Kl 
Limi  t=Corr 

ELSE 

Corr=3 

END  IF 
END  IF 


Plotting  of  the  data 


Kmax=K 

REDIM  D2(Knax+15) 

| 

!   positioning  the  data  on  the  plot  by  time  of  trace 
I 

Time=(Npoint2/Nplot  )*420 

Timedist=2  +  ((Npoint2+Npoint1  )/Nplot  )»420 

IF  Npoint1=0  THEN  Time=5 

IF  NpoinKE  THEN  Timedist=G 

End_pt=5cale_y*2G0 

D2(2  )  =  Kna;< 

D2(3)=INT(Temp  ) 

D2(4  >-INT<Time> 

D2(5  >-INT(Tinediat» 10) 
WINDOW  0,420  ,0,2G0 
GRAPHICS  ON 
l 

!  set  the  vertical  increment 
| 

Inc  =  ( (Num*Tminc*Rspd*End_pt )/Maxrng  ) 

D2(G)  =  INT(Inc*1000> 

I 
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4400 

4410 

4420 

4430 

4440 

4450 

44B0 

4470 

4480 

4490 

4500 

4510 

4520 

4530 

4540 

4550 

4560 

4570 

4530 

4590 

4G00 

4G10 

4620 

4630 

4640 

4650 

4660 

4670 

4680 

4G90 

4700 

4710 

4720 

4730 

4740 

4750 

47G0 

4770 

47C0 

4790 

4nrhr>> 
\JXJHJ 

4310 


4340 
4850 
4860 
4370 
4330 


!  compute  the  intensity  of  the  return,  move  to  the  proper 

I  coordinates  and  plot  the  appropriate  colored  block 

! 

FOR  K-1  TO  Kmax 

Z=Point(K)/Ilvl 

IF  Z>1  THEN  Z=1 

IF  Z<0  THEN  Z=0 

AREA  INTENSITY  Z.Z.Z 

Dis=INT(K*Inc  >+1 

MOVE  Timedist  ,Dia 

RECTANGLE  Tine  ,Inc .FILL 
NEXT  K 
! 

!   keep  an  account  of  the  trace  numbers  taken  on  the  plot 

! 

Nrec=Nrec+1 

Rec_num=INT( Nrec  5 

D2( 1  )=Rec_num 

Npoint 1 =Npoint 


save  routine 


if  function  key  is  set  then 
the  plot  will  be  saved 


IF  Cntrl-1  THEN 

ASSIGN  SFilel  TO  FilelS 
OUTPUT  OFilel  ,Rec_num5D2< » ) 

END  IF 


graphics  dump  of  plot  after  15  minute  intervals 


IF  Timedist>415  THEN 
Prt_dmp:    PRINTER  IS  701 
PRINT  " 
PRINT  " 
PRINT  " 

Plotnum=Plotnum! 1 
DUMP  GRAPHICS  $701 
FOR  1=1  TO  Kmax 

Ct(  I  )=Ct(  I  >»<(Temp+273>A2>"<  1  ./ 

Ct( I  )=Ct(  I  )»( 1 ./(A*G  )  ) 

Ct(I  )=Ct(  I  )*( 1  ./Pt  ) 

Ct(I  )=Ct(I)/Pl3_lng 

IF   CtC I  XI  .0E-r00   THEN   Ct(  I  )=1  .0E-1  00 
NEXT    I 
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4803 

4000 

!  CTA2  Computations  and  plots 

4910 

4920 

1 

4930 

GCLEAR 

4940 

VIEWPORT  15,120,15,80 

4950 

WINDOW  0,300,0,Exrng 

49G0 

AXES  30,Exrng/I5,0,0,30,Exrng/3 

4070 

CLIP  OFF 

4980 

CSIZE  2, .6 

4990 

L0RG  G 

5000 

FOR  M=0  TO  300  STEP  30 

5010 

MOVE  M,-Exrng/45 

5020 

LABEL  "1 .0E";((M/30)-10) 

5030 

NEXT  M 

5040 

MOVE  140,-15 

5050 

CSIZE  4 

50G0 

LABEL  "LOG  OF  CT"2" 

5070 

LORG  8 

5000 

FOR  M=0  TO  Exrng  STEP  INT<Exrng/3) 

5090 

MOVE  0,M 

5100 

LABEL  M 

5110 

NEXT  M 

5120 

LDIR  PI/2 

5130 

LORG  G 

5140 

MOVE  -40,Exrng/2 

5150 

CSIZE  4 

51G0 

LABEL  "RANGE  IN  METERS" 

5170 

LDIR  0 

5180 

LORG  4 

5190 

MOVE  150,Exrng+3 

5200 

LABEL  "TIME  AVERAGED  CTA2     ";Site$ 

5210 

CLIP  ON 

5220 

MOVE  300,0 

5230 

FOR  1=1  TO  Krcax 

5240 

Ypl  =  <  Maxrng/Kmax  )*I 

5250 

Pnt  =  (LOG(Ct(I  ))/2. 3025851  )+10 

52G0 

IF  Pnt>10  THEN  Pnt=t0 

5270 

IF  Pnt<0  THEN  Pnt=0 

5200 

Pnt=Pnt *30 

5200 

DRAW  Pnt  ,Ypl 

5300 

NEXT  I 

5310 

DUMP  GRAPHICS  $701 

-J  -'  -  u 

PRINT  " 

5330 

PRINTER  IS  CRT 

5340 

GCLEAR 

5350 

K  1  =  1  0 

53G0 

Lir<ut=Corr 

5370 

Noise=10*Corr 
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5380 
5390 
5400 
5410 
5420 
5430 
5440 
5450 
5460 
5470 
5480 
5430 
5500 
5510 
5520 
5530 
5540 
5550 
55S0 
5570 
5580 


IF  CntrI-1  THEN 

ASSIGN  ©Filel  TO  * 
END  IF 
Cntrl=0 


if  the   plot  is  to  be  saved,  the  file  is  created 
and  named 


IF  Save_plt=1  THEN 

CALL  File_init<Disc_address$ ,Nrec  , Filel $  ) 

Save_plt=0 

Cntrl=1 
END  IF 


start  the  next  IS  minute  plot 


GOTO  Again 


next  trace 


toggle  set  if  plot  is  to  be  saved 


5590  END  IF 

5800 

5810 

5520 

5830 

5640 

5650  GOTO  Sync 

5860 

5670 

5630 

5830 

5700 

5710  Svpt:  Save_plt»1 

5720       GOTO  Read_sig 

5730 

5740 

57S0   !  completion  routine 


5 -inn     n.  ,  ,  a  .  r-rnn 

5790  END 

58C0  ! 

5810  ! 

6000  !   SUBROUTINE  SECTION 

GS10  ! 

6020  ! 
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G030  SUB  Freq_init(Freq$,N$) 

G040   ! 

6050   ! 

G0S0   f     SETUP  OF  THE  HP3314A  FUNCTION  GENERATOR 

G073   ! 

G080   ! 

G090  *   Ap$-"AP" 

G100     Frq$="FR" 

SI  10     Nn$»"NM" 

G120     En$="ENM 

6130     Vo$-"VO" 

GU0     HzS-'KZ" 

6150     INPUT  "FREQUENCY  DESIRED  CkHZ)  (5  RECOMMENDED )" ,Freq$ 

S1G0     INPUT  "AMPLITUDE  DESIRED  (U  .  .1.5V  MAX)  "  ,Amp$ 

G170     IF  UAL(Awp$)>1  .5  THEN 

G180        Amp$-"1.5" 

6190        PRINT  "AMPLITUDE  OF  FUNCTION  GENERATOR  IS  1.5  V" 

6200     END  IF 

6210-     INPUT  "NUMBER  OF  CYCLES  PER  BURST  (INTEGER)  (100  RECOMMENDED  )"  ,N$ 

6220     OUTPUT  707; "M03" 

6230     OUTPUT  707; "SR2" 

6240     OUTPUT  707;Ap$&Amp$&Uo$&Frq$&Freq$&Hz$&Nm$&N$SEn$ 

6250  SUBEND 

6260   ! 

6270   ! 

6280   I 

6500  SUB  Init_ad200 

6510 

6520 

6530 

G540 

6550 

65S0     Ad_sel_code=17 

G573  Dunrny=READIC(  Ad_5el_code  ,3 

G580     URITEIO  Ad_5el_code,0;0 

65G0     CONTROL  Ad_sel_code ,0;  1 

GG00   SUBEND 

GG10   ! 

6S20   I 

6G30   ! 


INITIALIZATION  OF  THE  A-D  CONVERTER 
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7000 

7010 

7020 

7030 

7040 

7050 

70G0 

7070 

7080 

7090 

7100 

7110 

71Z0 

7130 

7140 

7150 

71G0 

7170 

7180 

7190 

7200 

7210 

7220 

7230 

7240 

7250 

72G0 

7270 

7280 

7293 

I       _  <_  O 

7310 


SUB  Set  time(Zone$) 


SET  THE  TIME  DATE  RECORDER 


PRINT  "WHAT  TIME  REFERENCE  ARE  YOU  USING?   INPUT: 
PRINT  "      1  FOR  UNIVERSAL  TIME" 
PRINT  "     2  FOR  LOCAL  TIME- 
PRINT  "     3  FOR  YOUR  OWN  CLASSIFICATION" 
INPUT  K 
IF  K-2  THEN 

ZoneS=»"( LOCAL)"  ' 
ELSE 

IF  K-3  THEN 
INPUT 

ELSE 

ZoneS 

END  IF 
END  IF 
IF  TIMEDATE<DATEC14  AUG  1984")  THEN 

INPUT  "ENTER  ""DD  MMM  YYYY" " " ,Date$ 

INPUT  "ENTER  " "HR:MIN:SC" " " ,Tine$ 

SET  TIMEDATE  DATE( DateS  )+TIME< TimeS  ) 

PRINT  DATE$<  TIMEDATE ) ,TIME$( TTMEDATE  ) 

Tstart=TIMEDATE 

T0=Tstart  MOD  8G400   ' 
END  IF 


'WHAT  IS  YOUR  TIME  REFERENCE"  , ZoneS 


(UTC) 


SUBEND 
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7500 

7510 

7520 

7530 

7540 

7550 

7560 

7570 

7530 

7590 

7G00 

7810 

7S20 

7G30 

7G40 

7G50 

7GG0 

7G70 

7G80 

7G90 

7700 

7710 

7720 

7730 

7740 

7750 

77G0 

7770 

7730 

7790 

7800 

7310 

7320 

7330 

i  o43 

73G0 

73G0 

7370 


SU8  Read_ad200( INTEGER  Dat(*>  BUFFER .Counts ,Tine$ ,0ff_ave ) 


INFOTEK  A-D  ROUTINE  SET  UP  FOR  EXTERNAL  TRIGGER 


Ad_sel_code=1 7 

! 

! INITIALIZATION  OF  THE  A-D  CONVERTER 
! 

OUTPUT  Ad_sel_code; "RESET" 

OUTPUT  Ad_sel_code; "INTERNAL" , "COUNT  "iCountS , "HOLDON' 

OUTPUT  Ad_5el_code;"DELAYON'\  "SELECT  Isl   end" ,Time$ 

OUTPUT  Ad_sel_code; "STATUS" 

ENTER  Ad_sel_code;Resp$ 

IF  Resp$=" "  THEN 

ASSIGN  @Ad200  TO  Ad_sel_code; WORD 


triggering  of  the  HP3341A 


i 


IF  Off_ave=0  THEN 

TRIGGER  707 
END  IF 

ASSIGN  @Buf  TO  BUFFER  Dat(»> 
TRANSFER  @Ad200  TO  @Buf;WAIT 
OUTPUT  Ad_sel_code;"  " 
OUTPUT  Ad_sel_code; "STATUS" 
ENTER  Ad_3el_code;Resp$ 

IF  RespSO" THEN  ' 

PRINT  "ERROR-  ";Re5p$ 
END  IF 
ELSE 

PRINT  "ERROR  DURING  INITIALIZATION  =";Re5p$ 
END  IF 
SUBEND 
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8000 

8010 

8020 

8030 

8040 

8350 

80S0 

3070 

8080 

8090 

8100 

8110 

8120 

8130 

8140 

8150 

81B0 

8170 

8180 

8180 

8200 

8210 

8220 

8230 

8240 

8250 

82G0 

8270 

3230 

8290 

8300 

8310 

8320 

8330 

8340 

8350 

33G0 

8370 

33G0 


3410 
8420 
8430 
3440 
8450 
84B0 
8470 
8430 


SUB  Plot_setup<  Nplot ,Site$ .Maxrng ,Scale_y ,  Exrng ,Zone$ ) 


!  SET-UP  OF  THE  TIME  PLOT  ON  THE  CRT 


Scale_y=Maxrng/Exrng 
GRAPHICS  ON 
LINE  TYPE  1 
VIEWPORT  15,120,15,80 
WINDOW  0, Nplot  ,0, Exrng 

AXES  Nplot/ 15, Exrng/ 15,0  ,0  ,Nplot/3,Exrng/3 
CLIP  OFF 
CSIZE  4..B 
LORG  G 

T1=TIMEDATE  MOD  86400 
Hrs=T1  DIV  3G00 
T2=T1  MOD  3G00 
Min-T2  DIV  G0 
Qtrhr-Min  DIV  15 
FOR  M-0  TO  Nplot  STEP  TNT(Nplot/3) 

MOVE  M,-Exrng/4S 

Qtrnin=Qtrhr»I5+<M»3/Nplot )*5 

IF  QtrNin=S0  THEN 
Qtrnin=0 
Hrs=Hrs+1 

END  IF 

LABEL  USING  "DO  , A  ,ZZ" ;Hra ; " : " ;Qtrnin 
NEXT  M 

MOVE  Nplot/2,-15 
LABEL  "TIME  "iZoneS 

LABEL  ORDINATE 

LORG  3 

FOR  M=0  TO  Exrng  STEP  INT(Exrng/3) 

MOVE  0,M 

LABEL  M 
NEXT  M 
LDIR  PI/2 
LORG  G 

MOVE  -Nplot/7,Exrng/2 
LABEL  "RANGE  IN  METERS" 

TITLE 

LDIR  0 
LORG  4 
MOVE  Nplot/2  ,Exrng+3 
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8490        LABEL  DATE$<  TIMEDATE  U  ""     ";SiteS' 

8500        CLIP  ON 

8510   SUBENO 

8520  ! 

8530  ! 

8540  ! 

9000  SUB  Attenuation(Freq ,Temp .Atten) 

9010  ! 

9020  ! 

9030  !  This  subprogram  calculates  the  attenuation  of  the 

9040  !  sound  in  air  based  upon  equations  in  Neff  1975 

9050  !  (source  of  subroutine:  Thesis  of  R.  Fuller) 

90G0  I 

9070  !   Variables 

9080  I 

9090  I  Atom_pres   -   input  atmospheric  pressure  in  mb 

9100  ! 

9110  !  Atten  -  attenuation  coefficient  of  acoustic  wave 

9120  ! 

9130  !  Att_max   -  Variable  in  program.   It  is  the  attenuation 

9140  !  at  the  frequency  of  the  maximum  attenuation 

9150  I  for  the  given  input  conditions. 

91G0  ! 

9170  !  F   -   Ratio  of  the  frequency  to  frequency  at  maximum 

9180  !       attenuation. 

9190  ! 

9200  !  Fmax   -  Frequency  of  the  maximum  attenuation 

9210  I 

9220  I  H  -  variable  used  In  the  integration  of  excess  attenuation 

9230  ! 

9240  !  Pstar   -  variable  used  in  intermediate  calculations 

9250  ! 

92G0  !  Tstar   -   variable  used  as  an  intermediate  in  calculation 

9270  I  of  attenuation. 

9280  ! 

9290  !  Uat_pres   -   Atmospheric  water  pressure  in  mb . 

9300  ! 

9310  INPUT  "Enter  the  atmospheric  pressure  in  mb" ,Atom_pres 

3323  INPUT  "Enter  the  atmospheric  water  pressure  in  mb" ,Wat_pre5 

9330  H=100*Wat_pres/Atom_pres 

9340  Tstar=( 1 .3*Temp+492  )/S1 3 

9350  Pstar=Atom_pres/1014 

9380  Fmax=( 1 0+GG00*H+44400*H»H  )*Pstar/Tstar " . 3 

9370  Att_max-.0078*Fmax*Tstar"(-2.5)»EXP(7.77*( 1-1/Tstar)) 

3380  F=Freq/Fmax 

9390  Atten-(Att_max/304.8)»((-.f8»Fr2+(2»F*F/(  1  +F*F  )  T2  )"  .5 

9400  Atten*(Atten+1 .74E-1 0*Freq*Freq )/4.35 

9410  SUBEND 
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9420  ! 

9430  ! 

9500  SUB  File_init(Di5c_address$,Nrec ,File1$) 

9510  1 

9520  ! 

9530  !   CREATE  THE  STORAGE  FILE  ON  THE  DISC  FOR 

9540  !   THE  REDUCED  DATA 

9550  I 

95B0  ! 

9570     INPUT  "ENTER  THE  REDUCED  DATA  OUTPUT  FILENAT1E  ",File$ 

9580     File!$=File$&Di3C_addres5$ 

9590     CREATE  BOAT  Filel $ , 1 80 ,400   • 

9B00     ASSIGN  ©Filel  TO  FilelS 

9610  ! 

9G20  SUBEND 
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APPENDIX  B 

SPEAKER  AND  ARRAY  ANALYSIS 
In  the  design  of  our  echosounder,  it  was  determined 
that  a  rapid  decay  time  was  required  to  obtain  accurate 
short  range  information.  Weight  restrictions  involved  with 
equipment  transportation  plus  maximum  response  at  high 
frequencies  led  to  our  decision  to  use  piezo  ceramic 
speakers.  The  Motorola  KSN  1005A  speaker  was  selected 
based  on  these  requirements  and  the  specifications  charted 
in  the  Motorola  catalog  [Ref.  16]  and  reproduced  below  in 
Figures  23  through  25. 

Nominal  Power,  Impedance,  and  Distortion  Ratings 
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Fig.  23.  Speaker  Ratings 
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Typical  Frequency  Response 
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Fig.  24.  Speaker  Frequency  Response 


Dimensions:  KSN  1005A,  KSN  1003A 


p  MMMfW  ^tf  * 


FOUR   218"  15  5mm|  DIAMETER  HOLES  EOUALLY 
SPACED  ON  A  3  94"  (100.1mmi  DIAMETER  BC. 


WEIGHT:  75  GRAMS 


Fig.  25.  Speaker  Dimensions 
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The  frequency  response  charts  indicate  that  a  maximum 
response  for  our  speakers  occurs  at  a  resonant  frequency  of 
5000  Hertz.  This  frequency  was  used  as  the  baseline  from 
which  all  our  measurements  are  made. 

Using  the  speakers  in  an  anechoic  chamber,  the  average 
e_1  voltage  decay  time  was  measured  to  be  approximately  900 
usee  (Fig.  26).  This  decay  time  translates  into  a  sound 
propagation  distance  of  just  over  15.0  centimeters  (at  STP) 
from  the  speakers.  Considering  our  requirements,  this 
speaker  is  ideally  suited  to  serve  our  purpose. 


At 


At=2.924  ms 
av=  25  mv 


/y\/\/VVV^^  •v./n.^-s^ 


T 


AV 


1 


Fig.  26.  Speaker  Decay  Time  Trace 
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B.   ACOUSTIC  ARRAY 

Our  next  consideration  was  the  echosounder  array 
pattern.  Ideally,  the  acoustic  sources  should  be  placed 
exactly  one  half  wavelength  apart.  At  a  frequency  of  5000 
Hertz,  this  would  require  spacings  of  3.4  centimeters  (at 
STP)  which  is  physically  impossible  for  the  speakers  we 
have  chosen.  The  closest  possible  spacing  is  7.62 
centimeters  between  sources  after  shaving  off  the  flange  of 
the  horn  (Fig.  25) . 

From  Kinsler,  et  al .  [Ref.  14],  the  equation  for  the 
directionality  factor  of  a  simple  line  array  is  derived  as: 


H(N,8) 


N 

sin 

( 



kd 

sin  9  ) 

1 

2 

N 

1 

sin 

( 

2 

kd 

sin  9  ) 

where 

k  is  the  wave  number  ( 27c/  a  ), 

d  is  the  distance  between  sources, 

N  is  the  number  of  sources,  and 

9  is   the  angle   measured  from   a  line  perpendicular  to 
the  array  to  the  direction  of  interest. 

However,  this  equation  assumes  simple  point  sources  which 
does  not  adequately  describe  our  speakers.  It  was  necessary 
to  couple  this  equation  to   the   directional   factor   for  a 
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piston  source   which  is   also  identified  in  Kinsler,  et  al 
[Ref.  14:  p.  108]  as: 


2   Ji  (  ka   sin  (  9  )  ) 
ka   sin  (  9  ) 


where 

k  is  still  the  wave  number, 

a  is  the  radius  of  the  piston  source, 

9  is  the  angle  measured  from  a  line  perpendicular  to 
the  array  to  the  direction  of  interest,  and 

Ji  (  x  )  is  a  Bessel  function  of  Order  1  with  an 
argument  x . 

These  two  equations  were  combined  together  to  form  the 
equation  of  directionality  for  a  linear  array  of  piston 
sources,  L(N,9),  by  simple  multiplication. 

L(N,9)   =   H(N,9)  •  D(9) 
This  equation  was   incorporated   into   a   computer  program, 
which  may   be  found   in  Reference   3,  and   used  to  generate 
linear  array  patterns  for  various  numbers  of  piston  sources 
(Figs.  27  through  31). 

It  was  concluded  from  these  plots  that  five  linear 
elements  would  provide  the  best  combination  of  forward 
directionality,  sidelobe  suppression,  physical  size  and  a 
relatively  low  cost.  Then,  in  order  to  enhance  both  array 
efficiency  and  symmetry,  we  settled  on  a  five  by  five 
element  array  design   with  7.26   centimeter  spacing  between 
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PRTTERN  FOR   3   ELEMENTS 
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Fig.  27.  Three  Element  Array    Fig.  28.  Four  Element  Array 


PRTTERN  FOR   5   ELEMENTS 


Fig.  29.  Five  Element  Array 


PRTTERN  FOR   B   ELEMENTS 


PRTTERN  FOR  7       ELEMENTS 


Fig.  30.  Six  Element  Array     Fig.  31.  Seven  Element  Array 
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speakers  in  both  the  vertical  and  horizontal  directions. 

After  verifying  the  manufacturers  polarity  designation 
for  35  speakers,  we  obtained  Lissajous  plots  for  each 
individual  speaker  (Fig.  32).  Based  on  the  speaker's  output 
to  input  voltage  ratios  as  illustrated  by  these  plots,  we 
were  able  to  rank  all  our  speakers  by  signal  efficiency. 
It  was  based  upon  this  criteria  that  we  selected  the  25 
most  efficient  speakers  for  the  array,  placing  the  best 
speakers  at  the  center  and  subsequent  ranking  speakers 
further  toward  the  sides  and  corners. 
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tea 

• 

spea! 

OUtDi 
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it     ^ 
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weaker 

Fig.  32.  Lissajous  Plot 
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The  25  selected  speakers  were  mounted  in  a  five  by  five 
planar  array  on  a  balsa  wood  insulated  bilayered  sheet 
metal  board.  After  wiring  all  the  speakers  in  parallel,  we 
surrounded  all  the  electrical  connections  and  speaker  backs 
with  two  3.0  centimeter  layers  of  foam  insulation 
sandwiching  a  1.0  millimeter  lead  sheet.  Then  the  entire 
array  mounting  was  enclosed  in  a  44  by  44  by  5  centimeter 
sheet  metal  box.  This  design  was  chosen  to  suppress 
virtually  all  acoustic  energy  propagating  out  the  rear 
hemisphere  of  the  array,  while  shielding  the  array  from  any 
external  electrical  interference  (Figs.  33  and  34). 
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Fig.    33.    Array  Photo 
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Fig.  34.  Array  Photo 


69 


APPENDIX  C 

ENCLOSURE  DESIGN 

Acoustic  echosounding  has  proven  to  be  an  extremely- 
useful  technique  for  probing  and  analyzing  the  lower 
atmosphere.  In  order  to  most  efficiently  utilize  the 
acoustic  waves  transmitted  and  later  received  by  this 
remote  sensing  method,  it  is  essential  to  have  an  efficient 
antenna  with  highly  directive  beams  and  strongly  suppressed 
sidelobes.  Antenna  design  becomes  increasingly  more 
important  in  a  noisy  environment  where  noise  pollution 
within  the  sidelobes  may  dominate  the  desired  signal  within 
the  main  lobe.  Hall  and  Wescott  [Ref .  12]  showed  that 
sidelobe  suppression  improved  with  higher  frequencies. 
Their  studies  showed  that  the  measured  90  degree  sidelobe 
suppression  ranged  from  38  dB  at  1  KHz  to  50  dB  at  5  KHz. 
Furthermore,  any  significant  improvement  in  sidelobe 
suppression  could  only  be  obtained  by  surrounding  the 
antenna  with  an  acoustic  energy  absorbing  cuff  or  shroud. 

In  an  effort   to   maximize   our   antenna  main   lobe  to 
sidelobe   power   ratio,   we   intend  to  operate  only  at  high 
frequencies   as    discussed    in    the    previous   section. 
Additionally,  we  have  designed  an  acoustic  energy  absorbing 
enclosure.     Many   designs   were    considered  based   upon 
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previous  research  in  the  field  of  echosounding  [Refs.  17 
through  20].  In  addition,  we  obtained  the  actual  acoustic 
beam  patterns  for  our  array  using  a  computer  program 
written  by  LCDR  Butler  [Ref .  21]  which  we  modified  for  our 
purposes.  This  modified  version  of  LCDR  Butler's  program 
may  be  found  in  Reference  3.  By  rotating  the  array  in  an 
anechoic  chamber,  we  were  able  to  produce  highly  accurate 
polar  plots  of  the  array  beam  patterns  (Figs.  35  and  36). 

5X5   ACOUSTIC   ARRAY   BEAM   PATTERN 


Jmike 


R-3.60  M         MIKE  VT=  8.2465       INPUT-5. 0  V         FREQ=5. 0  KHZ 
Fig.    35.    Polar  Plot   of    5   X   5   Array   Beam  Pattern 
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DIAGONAL  5X5  ACOUSTIC  ARRAY  BEAM  PATTERN 


R=3.60  M         MIKE  VT=  3.9497        INPUT=2. 5  V         FREQ=5. 0  KHZ 


Fig.    36.    Polar  Plot   of    5   X  5  Array  Beam  Pattern 

at   a   45°    Aspect 
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These  polar  plots  and  the  actual  coordinates 
corresponding  to  the  individual  data  points  confirmed  the 
computer  prediction  for  the  five  element  line  array  (Fig. 
29).  Additional  polar  plots  obtained  by  varying  both  the 
input  array  voltage  and  range  between  array  and  microphone 
further  support  our  claim  that  for  a  5  KHz  carrier 
frequency,  the  main  lobe  is  confined  to  a  divergence  angle 
of  20  degrees.  Since  it  is  our  aim  to  suppress  all 
sidelobes  and  utilize  solely  the  main  lobe,  we  chose  not  to 
taper  our  enclosure  as  most  previous  researchers  had. 
Rather  we  designed  the  enclosure  based  upon  the  dimensions 
of  the  array  itself  and  the  acoustic  beams  it  generated. 

Plywood  was  used  for  the  construction  of  the  enclosure 
and  provided  not  only  a  rigid,  inexpensive  framework,  but 
also  proved  to  greatly  attenuate  external  noise 
interference.  Anticipating  all  kinds  of  weather  conditions 
during  data  collection,  the  plywood  enclosure  was  first 
waterproofed  with  four  coats  of  marine  varnish.  Grooved 
joints,  caulking  and  weather  stripping  were  also  design 
considerations . 

A  millimeter  layer  of  lead  can  suppress  an  acoustic 
signal  as  much  as  40  dB  (Figs.  37  and  38).  About  7.0 
centimeters  of  corrugated,  egg-carton  design  foam  can 
suppress  a  signal  another  3  to  4  dB  (Figs.  37  and  39). 
Together  they  make  an  extremely  efficient  absorbing 
material  for  use  in  our  enclosure. 
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Fig.  37.  Response  Reference,  No  Insulation 
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Fig.  38.  Signal  Suppression  by  Lead 
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Fig.  39.  Signal  Suppression  by  Fo 


am 


Lead-lined  absorbing  foam  was  glued  to  all  inner 
surfaces  of  the  enclosure  with  two  1.0  millimeter  layers  of 
lead  overlapping  at  all  corners.  Strong  aluminum  brackets 
were  used  to  connect  the  four  side  panels  to  each  other  as 
well  as  to  the  enclosure  base  (Figs.  40  through  42). 
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Fig.  40.  Enclosure  Photo,  Fully  Assembled 
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Fig.  41.  Enclosure  Photo,  Interior 
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Fig.  42.  Enclosure  Photo,  Base  with  Array 
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APPENDIX  D 

ECHOSOUNDER  OUTPUT 

Sixteen  echosounder  output  traces  were  included  to 
exhibit  typical  atmospheric  activity.  Many  of  the  plots, 
such  as  the  14:45,  15:30,  and  the  17:30  of  26  April  and  the 
10:00,  11:15,  13:00,  14:15,  15:15,  16:15,  and  the  17:15  of 
27  April,  have  convective  plumes  which  are  prevalent 
whenever  a  heat  flux  between  the  surface  and  atmosphere 
exists.  The  plots  of  18:15  and  19:15  on  the  26th  of  April 
and  of  18:00,  18:45  and  20:00  on  the  27th  of  April  clearly 
show  the  passage  of  the  neutral  event  which  is  encountered 
when  the  atmospheric  and  surface  temperature  difference 
becomes  negligible.  Finally,  the  plot  of  16:30  on  26  April 
can  be  associated  with  strong  winds  which  exhibit  a 
somewhat  uniform  return  for  all  altitudes  across  the  entire 
15  minute  interval  of  the  trace. 

The  horizontal  axis  labels  of  the  Ct 2  plots  are  small 
and  may  be  hard  to  read.  This  axis  is  a  log  scale 
beginning  with  l#10_io  at  the  left  side  of  the  plot  and 
ending  with  1»10°  at  the  right  side  of  the  plot.  Each  tick 
mark  moving  left  to  right  along  this  axis  represents  an 
integer  increase  in  the  exponent  of  10. 
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